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Abstract. In hybrid inflationary models, inflation ends by a sudden instability 
associated with a steep ridge in the potential. Here we argue that this feature can 
generate a large contribution to the curvature perturbation on observable scales. This 
contribution is almost scale-invariant but highly non-Gaussian. The degree of non- 
Gaussianity can exceed current observational bounds, unless the inflationary scale is 
extremely low or the hybrid potential contains very large coupling constants. Non- 
linear effects on small scales may quench the non-Gaussian signal, and while we find no 
compelling evidence that this occurs, full lattice simulations arc required to definitively 
address this issue. 

Note added: We now believe that nonlinear effects will invalidate the original 
computation in this paper essentially instantaneously after the short-wavelength 
modes reach the minimum of their potential. This means that the mechanism 
described in this paper will not lead to appreciatable curvature perturbations on long 
wavelengths, and no useful constraints on hybrid inflation will result. We have 
inserted a brief calculation on p2 of this manuscript to explain this fact, but have 
otherwise left the manuscript unchanged. 



Keywords: Inflation, cosmological perturbation theory, physics of the early 
universe, quantum field theory in curved spacctimc. 
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Version 2 of this paper withdraws the claim of large / NL 

After further review, we have come to the conclusion that nonlinear effects cut off 
the growth of long-wavelength modes nearly instantaneously after the short- wavelength 
modes reach the minimum of their potential. This means that even the small difference 
in timescale explored in Section 5 is not sufficient to guarantee the growth of appreciable 
long-wavelength fluctuations. Hence we believe that the mechanism originally described 
in our paper does not lead to the growth of curvature perturbations on large scales. 
Below, we present an argument which illustrates this fact. 

We isolate the waterfall field x and assume it has a potential of the form: 

V(x) = ^(x 2 -v*)=V -^ + ^ (1) 

where fi 2 = Xv 2 . (Note \x has a different meaning in this note than in the rest of the 
paper). The equation of motion for \ is 

X - V 2 X - A + Ax 3 = 0, (2) 

and holds pointwise. We are really interested in the equation of motion for a related 
quantity, namely x "smoothed" on some length scale L. Let us denote this "smoothed" 
version of x( x ) by S£,[x](x). There are a few different ways to define this smoothed 
variable, such as convolution with a transfer function or via a cutoff in Fourier space, 
but all we really need to proceed with the argument is that the smoothing is linear, so 

sl[c\Xi + C 2 X2]0) = CiS L [xi](x) + c 2 s L [x2](x) (3) 

for all constants C\ and C2, and that the smoothing is a projection^], so 

sl[sl[x}] = s l [ X } (4) 

Since the smoothing process is a projection, we can use it to decompose x m to a 
"smoothed," long- wavelength piece xl, and a short- wavelength piece xs as follows 

X = Xl + Xs (5) 

where xl(x) = s L [x](x) and xs = X ~ Xl- 

We can now study how the short-wavelength dynamics affects the long-wavelength 
ones. Applying the smoothing process sl to the equation of motion (jSJ), one finds 

XI - y?XL + Ax| + AA = (6) 

where 

A = s L [ X s ]-xl (7) 

We have dropped the gradient term, since L is much larger than the other scales in 
the problem. Equation (jBJ) means that the smoothed field x has the same equation of 

| This property holds only approximately for convolutions, but the approximation will be sufficently 
accurate for large L. 
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motion as the unsmoothed field, up to an additional effective term A A. To estimate this 
term, we use the decomposition (jSJ) and the properties ([3]) and (jl]) to find 

A = 3 X 2 L s L [xs] + ZXls l [x 2 s\ + «x[x|] (8) 

The first term must vanish because of the definition of the decomposition (jSJ). The 
last term may be nonzero, but it seems likely that it should also vanish because of the 
X — > ~X symmetry in the problem. This leaves the middle term. Once the short- 
wavelength modes have reached the minimum of the potential, Xs = 1,2 Hence, once the 
short scale- wavelength modes have reached the minimum we have 

A = 3v 2 XL (9) 

Plugging this in to ([6]) results in 

Xl + 2/i 2 XL + Axi = (10) 

The short-wavelengthe modes have caused the appearance of an effective mass term in 
the equation of motion for xl- This term cancels the tachyonic mass term present in 
the action, and stabilizes the long-wavelength modes. This means that by the time the 
short-wavelength modes have reached the minimum of the potential, the long- wavelength 
modes should cease evolving essentially instantaneously. 

It is perhaps interesting that if only quadratic terms were present in the action (1), 
then the equations of motion for x would be linear and the tachyonic mass term would 
be preserved in the xl equation of motion. However, for inflation to end the potential 
must have a minimum, and the nonlienarities required to introduce the minimum also 
instantaneously cut off the growth of long-wavelength modes once the short-wavelength 
modes reach the minimum of the potential. 

The correct calculation 

The SN formalism can still be used to calculate the contribution to ( from the hybrid 
transition, once we account correctly for the short scales modes. When we smooth the 
universe on large scales, we must remember that although the waterfall field x averaged 
on these scales is extremely small, the kinetic energy density of the x field averaged on 
these scales has a contribution from modes of all scales, and indeed is dominated by the 
shortest scale modes present. The time taken for the universe smoothed on large scales 
to transition from a flat initial hypersurface to a kinetically dominated comoving final 
hyperspace is therefore not given by equation (46), as we had calculated, but equation 
(46) with x*, the initial field smoothed on large scales, replaced by <r*, the RMS value of 
the initial x field, with all modes included. This quantity is dominated by the shortest 
scale modes present in the problem, and can be read off from equation (43) using (44) 
with k given by fc s h or t- From the point of view of the large scale modes a* is almost a 
constant, but has a small cosmic variance given by x*- The induced variance in N and 
therefore (, due to this cosmic variance in <r* is therefore given by the expression 

5N = ( = -^L X *, (11) 



Non-Gaussianity constrains hybrid inflation 



4 



where N is given by equation (46) multiplied by H and with x* replaced by cr* 
as discussed above. Therefore, the curvature perturbation on a particular scale 
is proportional to the field perturbation averaged on that scale, and is therefore 
exponentially suppressed on the largest scales. 
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1. Introduction 

Recently, a rapid accumulation of data has made clear that the orthodox scenario of 
structure formation — in which small fluctuations oscillate under the influence of gravity 
within a almost-uniform primeval plasma — is in excellent agreement with observation. 
The initial conditions for these fluctuations are determined by the primordial curvature 
perturbation, £, usually supposed to have been synthesized during an earlier cosmic 
epoch, perhaps during inflation or ekpyrosis. Many versions of inflation are characterized 
by slow, smooth evolution and yield very Gaussian fluctuations. In other versions 
evolution is a dramatic process, allowing larger non-Gaussianities. 

In this paper, we show that significant power and non-Gaussianity can be generated 
when a set of inflationary trajectories fall from a ridge in field space. Such a process 
is a crucial element of the hybrid inflationary scenario, introduced by Linde, in which 
inflation ends through a sudden instability [H |2l [3]. Quantum fluctuations lead to slight 
differences in field values at the top of the ridge, causing different regions of the universe 
to follow one of a narrow bundle of trajectories. Trajectories near the edge of the bundle 
are ejected from the ridge early in their evolution, whereas those in the core remain on 
the arete longer. The dispersion generated in this way induces a variation in expansion 
history ("5iV") from trajectory to trajectory within the bundle, resulting in a nearly 
scale-invariant but highly non-Gaussian spectrum of density perturbations. A careful 
analysis of this process reveals the surprising fact that the final dispersion in the bundle 
depends only weakly on the initial dispersion. Hence falling from the ridge can have 
a large effect on the final curvature perturbation, even if the bundle of trajectories is 
initially very tightly focused. 

Ultimately these contributions to the power spectrum can be traced to the 
conversion of isocurvature modes into adiabatic ones during descent from the ridge. 
In the early days of the inflationary paradigm, it was assumed that the curvature 
perturbation was primarily determined by quantum fluctuations in the inflaton field 
"freezing out" as successive physical scales left the cosmological horizon. Later, when 
inflationary models containing two or more fields were constructed, it was understood 
that the curvature perturbation could grow or decay as neighbouring horizon volumes 
diverged along adjacent but inequivalent field-space trajectories. The evolution is caused 
by light isocurvature degrees of freedom which play no dynamical role while cosmic 
microwave background (CMB) scales are leaving the horizon, but later come to influence 
the expansion history of the universe. Although models with isocurvature effects are 
more complicated, they admit a richer phenomenology. In models with canonically 
normalized fields the reprocessing of isocurvature modes can lead to a significant 
non-Gaussian component of £ [U |5j [6]. In the present case, the isocurvature modes 
are massive while CMB scales leave the horizon and are suppressed by cosmological 
expansion. Nevertheless, falling from the ridge enormously amplifies these perturbations 
and can produce a significant effect. 

Our analysis applies to the original hybrid scenario and many related models. The 
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qualitative conclusions may apply during a phase of two-field ekpyrosis [7J El EJ QUI 
[Til H~2] . which also depends on dispersion near a ridge. The inflationary production 
of non-Gaussianity in dispersing models has been source of interest in recent years. 
Alabidi determined the non-Gaussian yield in a hybrid model, beginning from a range 
of locations near but not coincident with the ridge [13]. Byrnes et al. performed a 
similar analysis, showing that the initial conditions could be tuned to produce a large 
effect [TU [15]. Related studies have appeared in the literature [HI \T7\ . 

The generation and evolution of fluctuations in hybrid inflation was studied by 
Randall, Soljacic & Guth [18], and later Garci'a-Bellido, Linde & Wands [HI [20]. More 
detailed calculations were undertaken by Copeland, Pascoli & Rajantie [21]. The field 
whose instability ultimately ends inflation is termed the waterfall field and at early 
times its potential is designed so that it is very heavy, causing its perturbations to 
decay. Only fluctuations orthogonal to the waterfall make an appreciable contribution 
to the number of e-folds, N, required to reach a subsequent surface of uniform energy 
density. The result is a spectrum of density perturbations which can be evaluated 
by ignoring the waterfall field. In the simplest scenario, where only one canonically 
normalized field <fi is relevant, the single-field formula 5N ~ H8(j>/\(j>\ applies. There is an 
added complication, because the waterfall field rolls very slowly at the epoch of ejection. 
Therefore fluctuations in N may be large on scales leaving the horizon at that time, and 
one must verify that there is not an unacceptable synthesis of topological defects or 
primordial black holes [20J. The original hybrid model produced a blue spectrum of 
perturbations, but a larger class of so-called "hill-top" models [221 E3] yield red spectra 
compatible with CMB constraints [24]. 

In this paper, we point out that the traditional calculation of curvature 
perturbations in this model should be augmented with a new contribution sourced by 
the waterfall transition. While the trajectories remain on the ridge, fluctuations in 
5N are well-described by conventional perturbation theory applied to the transverse 
directions, as explained above. At the epoch of ejection, however, fluctuations in the 
waterfall field become important, no matter how small they are. To see that this can 
be so, note that 8N measured between the central trajectory (which does not leave 
the ridge) and any other trajectory (which falls into a non-inflating minimum) tends to 
infinity. Thus, even exponentially suppressed fluctuations can be amplified. It follows 
that a tiny fluctuation between neighbouring horizon volumes may ultimately give rise 
to a large density fluctuation. 

The basic mechanism by which large non-Gaussianity is generated can be illustrated 
using an inverted simple harmonic oscillator. Suppose we consider a particle of unit mass 
and position x, for which the Lagrangian 

L = \ ±2 + (12) 
describes an inverted harmonic oscillator of natural frequency u>. We wish to compute 
the time t required for the particle to roll down to some specified final position, xp, as 
a function of its initial position x . We assume that lox? 1, xq <C Xp, and that the 
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particle has zero initial velocity. At late times x(t) ~ (x /2)e UJt , so 

t(x ) = -\np£ (13) 

CO \Xo\ 

Next we study the statistics of t for an ensemble of particles in which x is Gaussian 
distributed with variance a 2 , assuming cocr <C 1. The mean rolling time is 



^2 . . . . . 

'Xp 



roo e -x Q /2a l o x 7 + 31n2 / 

(t) = / — ^= In - — - dx = - h In 

J-oo V27TO- to \x \ to 2 V 



a 



(14) 



where 7 is the Euler-Mascheroni constant. The average time (t) depends only weakly 
on the initial distribution width, a, with narrower distributions taking longer to roll off 
the arete. More surprising results obtain if we compute higher statistics of t, such as its 
variance 

°l = ((t ~ (t)) 2 ) = ^ (15) 



third moment 



/i 3 = ((t-(t)) 3 ) = ^ (16) 
and fourth moment 

7vr 4 

^ = «* -<*»*> = 647 (17) 
where ((z) is Euler's zeta function and £(3) ~ 1.202 is Apery's constant. Unlike (t), 
these moments are completely independent of the initial variance a 2 , and depend only on 
the curvature of the potential. Therefore, no matter how small the variance associated 
with xq, the variance in arrival times is the same. Moreover, one typically measures 
departures from Gaussianity by computing quantities such as the skew and 
excess kurtosis n^/af — 3. In our case these are pure numbers, of order unity, which are 
completely independent of 00. Therefore, our purely Gaussian initial distribution always 
becomes highly non-Gaussian through the process of rolling off the hill. 

This toy example can be translated into the inflationary context using a simple 
dictionary. Roughly speaking, the coordinate x represents the waterfall field which 
mediates the end of hybrid inflation, the end point Xp represents the location of the 
reheating minimum, and the time t represents the number of e-foldings required to 
reach reheating, at which point inflation ends. The initial Gaussian distribution of 
release points models the quantum flucutations in the waterfall field just before the 
hybrid transition. Fluctuations in the e-folding history are related to the curvature 
perturbation, (. Therefore, the moments of t represent the moments of (. The fact 
that they are large suggests that large non-Gaussianity can be generated by the hybrid 
waterfall. Of course, the full inflationary story is much more complicated, and we will 
support our conclusions with semi-analytic arguments and numerical calculations using 
realistic forms of the inflationary potential. However, the basic mechanism is sufficiently 
general that its essence is captured by the inverted harmonic oscillator. 

This mechanism generates power on all scales. In the traditional picture of hybrid 
inflation, effects occuring near the waterfall transition could only generate significant 
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power on short scales. For example, Randall et al. and Garcfa-Bellido et al. identified 
large fluctuations which collapse to black holes or topological defects. These are 
generated on the horizon scale at the epoch of ejection, yielding a spike in the power 
spectrum at small wavelengths. In contrast, a careful analysis of the process described 
above indicates that its power spectrum is nearly scale-invariant. This is possible 
because the final variance between expansion histories depends very weakly on the initial 
dispersion of release points, as illustrated by the toy model above. During inflation, long- 
wavelength fluctuations in the waterfall field are inevitably produced, although these 
are strongly suppressed by cosmic expansion. In the picture we employ, these tiny 
fluctuations are encoded in the width of the bundle of field-space trajectories followed 
by cosmological-scale patches of the universe. After the hybrid transition, the tilt of 
the power spectrum can be determined by comparing the dispersion of the bundle of 
trajectories followed by large-scale patches with that of horizon-scale patches. Since the 
final dispersion of trajectories is very weakly dependent on the initial dispersion, the final 
dispersion in each bundle is roughly the same. This is the hallmark of a scale-invariant 
process. 

Could this new component in the curvature perturbation become observationally 
detectable? If the growth of perturbations during the waterfall is not curtailed it will 
be amplified indefinitely, yielding a bispectrum in conflict with observation. However, 
any hybrid model automatically contains a cut-off on the amplification which can be 
achieved. While traversing the ridge, small fluctuations push horizon-scale patches to 
one side or the other, which roll down and reheat into radiation. This so-called tachyonic 
preheating [251 [26J [27] is non-perturbative, and does not admit a description in terms of 
a homogeneous coarse-grained background field with small fluctuations. Beyond some 
limiting time, of order that required for the field to reach the minimum, the evolution of 
large-scale modes becomes dominated by radiation produced during the reheating phase 
rather than the inflationary potential. Because the ridge in field space no longer drives 
trajectories to diverge, the growth of isocurvature modes from the waterfall ceases. 

Ultimately, numerical simulations may be required to determine whether a 
significant effect can be generated before perturbation theory breaks down. One might 
have expected that a very long time would be required before isocurvature fluctuations 
from the waterfall could compete with the adiabatic mode. In fact, our analysis suggests 
that the field smoothed on large scales typically reaches the reheating minimum only 
a fraction of an e-fold after the field smoothed on the horizon scale. Whether a given 
model can generate acceptably Gaussian fluctuations depends on a subtle competition 
between the perturbative growth of isocurvature modes and a non-perturbative cut-off 
from tachyonic preheating. In our opinion, it is not clear that the outcome can be 
decided on the basis of perturbative calculations. Our aim is to show that there is 
serious interest in obtaining the answer. We hope that the quantitative intricacies of 
the competition can eventually be settled using lattice simulations. 

The structure of this paper is as follows. In £J21 we fix our notation by reviewing 
familiar material concerning inflationary perturbations and the SN formalism. In §3] we 
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give a phenomenological description of evolving field fluctuations near a sharp ridge, 
and use this to derive some approximate expressions for the spectrum and bispectrum 
sourced in the curvature perturbation. In £JH we discuss the validity of our calculation. 
We work in natural units where c = h = 1, and the Planck mass associated with 
Newton's gravitational constant is set to unity, Mp = (87rG) -1 / 2 = 1, although Mp will 
occasionally be restored to clarify the relative magnitude of terms. 

2. The 5N formalism 

Experience has shown that an especially convenient accounting of isocurvature modes 
can be obtained using the Ll 5N formalism," introduced by Starobinsky [28] and Stewart 
& Sasaki [29], and extended beyond linear order by Lyth & Rodriguez [5]. The 5N 
formula is an example of the "separate universe picture," according to which distant 
Hubble volumes evolve like separate unperturbed universes. 

The comoving curvature perturbation, (, measures fluctuations in the expansion 
history between spatially disjoint regions of the universe. According to the separate 
universe picture it can be evaluated using the 5N formula, 



where t* labels an initial spatially flat hypersurface, and t c labels a subsequent uniform 
energy-density hypersurface on which the energy takes a prescribed value p c . 

Eq. (JTBl is the lowest term of an expansion in powers of the dimensionless 
gradient k/aH but is otherwise non-perturbative. Analytic calculations can frequently 
be simplified by constructing a Taylor expansion of iV(p c ,t*,x) around an arbitrary 
fiducial trajectory. When inserted in expectation values, the high order terms in 
this expansion generate contributions which grow with volume, describing a delicately 
shifting pattern of correlations in the large- volume limit. In analogy with the growing 
ultra-violet ("loop") contributions which govern correlations in quantum field theory 
as one samples fluctuations of increasing energy, these terms have been interpreted 
as "classical loops" or "c- loops" [301 EI], an d it has become popular to frame 5N 
calculations in this language. In this paper we do not make use of a loop expansion, 
but apply the non-perturbative definition, Eq. ({TBI , directly. We believe this choice has 
important advantages. First, issues associated with convergence are avoided. Second, 
the calculation need not be artificially truncated at a low order in the expansion. Third, 
it is unnecessary to introduce an infra-red regulator or "factorization" scale (without 
physical significance) to make intermediate steps of the calculation finite. In quantum 
field theory the loop expansion is unavoidable owing to our ignorance of ultra-violet 
physics. In contrast, in many models the infra-red physics associated with Eq. (TIB"]) is 
under reasonable control and in these cases we believe the loop expansion is, at best, 
an avoidable complication. In unfavourable cases it may be rather misleading. 

In the following sections, we supply expressions which allow N to be calculated in a 
hybrid model. Eq. f fl8|) relates the statistical properties of N to those of the observable 
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quantity, £, by the rules 

(CO = (N 2 ) - (N) 2 , (19) 
(CCC) = (N 3 ) - 3(N) (N 2 ) + 2(N) 3 . (20) 

Similar formulae can be written for higher-order correlation functions. 

Eqs. (fT9l- (l20l give the variance and skew of a collection of independent spacetime 
volumes smoothed on some characteristic scale, L. They are not spatially dependent. 
To obtain the real- and Fourier-space correlation functions which can be compared with 
observation, one makes a second use of the separate universe picture. This implies that 
any two spatially disconnected volumes will be uncorrelated. Therefore, 

(U y ) = a 2 (L)L 3 5(x-y) } (21) 

where we have written (CC)z = ° 2 In Fourier space, the power spectrum P{k) is 
defined by 

(Ck 1 Ck 2 )=P(^i)(27r) 3 5(k 1 + k 2 ), (22) 

and after making the identification k ~ L^ 1 we can conclude that P(k) = a 2 (k)/k 3 . 
The bispectrum, -B(ki, k 2 , ka), can be obtained by similar means. It satisfies 

(C kl Ck 2 Ck 3 ) =5(ki,k 2 ,k 3 )(27r) 3 (5(k 1 + k 2 + k3). (23) 

It is conventional to use momentum conservation to make B a function of the magnitudes 
{k±, k 2 , k 3 } alone. We write (CCC)l = oc(L). Then, 
/> > > v oi(L) 

(CxCyCz) = -y- 

x {L 3 L 3 5(x - y)5(x - z) + L 3 Ll6(y - x)<J(y - z) + L 3 L 3 5(z - x)5(z - y)} . (24) 
Comparing Eqs. (I2"3"j) and ffM§ . it follows that the bispectrum can be written 

B - % ftf' (25) 

which is the 'local' form typically generated by superhorizon evolution. This is a natural 
consequence of Eqs. (l2Tj) - (l2~4"j) . which require that correlations are local in real space. 

To measure the amplitude of three-point correlations it is conventional to use the 
/nl parameter, defined by 

B(h, k 2 , h) = p NL {p(h)P(k 2 ) + P(h)P(k 3 ) + P{k 2 )P(k 3 )}. (26) 



Eqs. (1251) and (1261) imply that /nl satisfies 



5 a(L) 5 (CCC) 



/nl - — ., T , - — tttto • (27) 



18a 4 (L) 18 (CC) 



This formula can be derived from many alternative constructions (see Ref. |32j). 
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3. Conversion of isocurvature fluctuations at the waterfall 

Consider any region of field space described by a collection of light scalars (f) a and a 
waterfall field Xi evolving according to the action 

S = -l [ d A x^\drd ( j> a + (d X ) 2 + 2V(r,x)}, (28) 



2 . 

where (dx) 2 = d a xd a X an d lower-case Latin labels {a, b, . . .} index spacetime 
dimensions. The potential V(<p a ,x) should be suitable for realising hybrid inflation 
but is otherwise arbitrary. 

Near the point at which the waterfall field, Xi becomes tachyonic, many hybrid 
potentials can be well-approximated by 

V = V(<p,...)-\<p X 2 + .... (29) 

where V involves only fields orthogonal to x- The field (p is a single direction in field 
space which behaves as an order parameter controlling the onset of the waterfall, and 
A is a coupling with the dimensions of mass. If the density perturbation is dominated 
by adiabatic fluctuations among these fields, then the spectral tilt, n s , depends on V 
alone. It is given by 

n s - 1 = -6e + 2r], (30) 

where the slow-roll parameters are defined by 

M ^tX and ^Mlt. (31) 
V V 



A prime ' denotes a derivative with respect to the adiabatic direction (p. Linde's original 
formulation of hybrid inflation yielded a blue spectrum [1] , which is now known to be 
incompatible with observation. In hybrid models e is typically negligible, so acceptable 
V generally require r\ to be slightly negative, satisfying r\ « —0.02. Such models 
are sometimes described as "hill-top" potentials [23]. Potentials with a negative rj 
complicate the higher-order terms necessary to guarantee a consistent model, but are 
conceptually no different to the original hybrid proposal. We will use an explicit example 
of such a potential in §3.21 



3.1. Analytic approximation 

The dynamics of Eq. fl29l) can be complicated. Our analysis relies on certain simplifying 
approximations, and closely parallels that of Copeland et al. [2T]. We ignore Hubble 
damping and assume that <p is practically constant. These will be acceptable during 
ejection from the ridge, provided that ip is rolling slowly and the trajectory is ejected 
over a timescale much shorter than a Hubble time. Ultimately, we will justify these 
approximations by checking our results numerically. For this purpose we use both the 
effective potential, Eq. ( T2"91 . and global completions which describe the descent into a 
reheating minimum, to be discussed below. 
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It is convenient to work in Fourier space, where \ can be decomposed into a series of 
coupled oscillators Xk- We may freely choose coordinates so that the effective potential 
for x becomes tachyonic at if = 0, and set t = at that time. For t > the evolution 
of each % k is governed by a time-dependent equation of motion, 

Xk + (A: 2 -^ 3 t)xk = 0, (32) 

where the linear growth with t is a consequence of our assumption that <fi is practically 
constant. The overdot indicates a derivative with respect to t, and \x has dimensions of 
mass, 

// = 2A0| t=o . (33) 

The general solution is 

Xk(*) = A k Ai{fit - k 2 /^ 2 ) + 5 k Bi(//t - k 2 /fx 2 ), (34) 

where and B k are arbitrary coefficients, and Ai(x) and Bi(x) are Airy functions of 
the first and second kind, respectively. For x < these functions are oscillatory, whereas 
for x > the first Airy function decays and Bi(x) grows exponentially. This growth 
corresponds to an exponential rise in occupation number of the low-lying modes of X- 
It describes a rapid reordering of the spatial field configuration, during which gradients 
may become large. We will return to this issue in £jH 

For a mode of wavenumber k, exponential growth begins at the time tk = k 2 / /i 3 . It 
follows that the phase transition proceeds by spinodal decomposition, beginning with the 
zero mode and extending to higher A;-modes at successively later times. The coefficients 
Ak and B k can be expressed in terms of Xk(^fe) and Xk(^fe), giving 

3 2/3 3I/3 

A k = __r(2/3)xk(tfc) - — r(l/3) Xk (4) (35) 

B k = ^!r(2/3) Xk (4) + ^r(l/3)*k(t*)- (36) 

Our interest lies in a description of the process by which spatially disjoint spacetime 
regions depart from the arete. It is at this point that our analysis departs from that 
of Copeland et al, whose interest lay with the smallest possible scales. Consider an 
ensemble of spacetime volumes passing over the waterfall, smoothed on a lengthscale of 
order k~ x . The mutual scatter in field value between the members of this ensemble can 
be computed, for which it is a good approximation to include only the growing mode. 
We find 

a 2 x (k, t) = -L I p B { k ') Bi[/z(t - t k ,)] 2 k' 2 dfc', (37) 
^ Jo 

where we have filtered out modes of wavelength shorter than 1/k, and introduced the 
power spectrum of Bk, 

(B* k B k ) = P B (k)(2n) 3 5(k-k'). (38) 
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Eq. ( 13 7p is valid when the relevant fc-modes are already rolling, and hence applies 
provided k < ^/ /i 3 t. We have set an implicit infrared cutoff scale to zero. It transpires 
that Pb(/c) is a very blue function of k, so this is typically an excellent approximation. 

Determining Pb(/c) requires information about the quantum fluctuations in xljJ 
Immediately prior to the tachyonic transition, \ carries a power spectrum imprinted by 
the preceeding epoch, during which its mass was large. It is a reasonable approximation 
to determine the fluctuations synthesized during this era by setting each k-mode in a 
Bunch-Davies vacuum state, with large constant mass M. Well inside the horizon, each 
mode is normalized to the corresponding Minkowksi space oscillator. Once k/a < M 
the mass term will dominate its evolution, and each mode decays as Xk oc a -3 / 2 owing 
to Hubble damping. For super- horizon modes, where k/a ^ H, this leads tdj] 

«^=2Mf) 3 < 2 ^< k - k '> 

M / f-f \ 3 

(x k X k '> = y (j-J (27r) 3 <5(k-k'), (40) 

where fc* is the comoving wavenumber of the mode leaving the horizon at the time of 
evaluation. These correlators scale like a -3 , because the wavenumber k ± scales like a as 
inflation proceeds. Employing fl35|) and ( I4TJ]) yield^jj] 

2" 



P B {k) U 



8Mfc3 



3 l/3 r (2/3)2 + 3-l/3 r(1/3) 2 



where C aggregates the numerical constants. It is at least of order unity, and can be 
much larger if /i <C M. 

For large x, the growing Airy function Bi(x) takes the asymptotic form 

Bi(x) ~ -^=x- l/i exp ( -x 3/2 ) . (42) 



hr r V3 

At late times, this implies that the scatter among field values is well-approximated by 



ch 3 r k i /4r . ,i 3 / 2 



f In their analysis, Copeland et al. noted that modes are light at the tachyonic transition and assumed 
that each mode was in a massless Bunch-Davis vacuum state at time tk- In this state, the correlation 
functions satisfy 

<XkXk') = ^(27r) 3 5(k-k') 

(x k Xk'> = ^(27r) 3 5(k-k'). (39) 

This would be a good choice whenever the effective mass M of \ varies adiabatically, or whenever 
M /M 2 <C 1. Very close to the tachyonic transition, M 2 ~ t and evolution will not be adiabatic. 
| See | Appendix A| for details. 

§ We have evaluated the correlators (Xk(ffc)Xk'(ifc)) and (Xk(^)Xk'(^)) a t a common time, t = Q, 
rather than tk ■ The error we commit in this approximation is negligible, since the e- foldings of expansion 
between t = and t — tk is of order (k / k*) 2 (H j fi) 3 <C 1. 
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When t > ^ the fc-dependence becomes trivial: all modes behave collectively, evolving 
coherently with the background field. Moreover, the fc-integral can be performed at 
once. In this limit, we find 

<^ = §^iwM^ t]3 ' 2 ) ■ <44) 

Eq. (|44p implies that the scatter within an ensemble of spacetime volumes grows 
according to the equation of motion for the background field in a single volume, 

X* / 2 r ^3/2 



^w^ ex H3 M/ J- (45) 

This justifies our application of the separate universe picture, but the preceding Fourier 
space analysis is required to determine the effective release point, x*i as a function of 
scale. This can be obtained from inspection of Eq. ( |44p . 

3.2. The statistics of N 

Once the waterfall potential has become tachyonic, % rons down the ridge. Its 
velocity increases until the universe becomes dominated by the kinetic energy x 2 /2. 
When this phase of "kinetic domination" is achieved, comoving hypersurfaces are 
practically determined by hypersurfaces of constant kinetic energy, and ( ceases to 
evolve. According to Eq. f l4"5j) . the distribution of initial values, x*i will l ea d to a 
spread in arrival times at the kinetically dominated epoch. We write the perturbation 
in expansion history associated with this spread as ( w . It is this fluctuation which we 
propose must be included when calculating the properties of perturbations generated in 
a typical hybrid model. The following perturbative analysis will be trustworthy only if 
kinetic domination is reached. In §H] we discuss the possibility that a non-perturbative 
effect, tachyonic preheating, can quench the dispersion of trajectories before kinetic 
domination is attained. 

The transit time to kinetic domination, t(x*), can be estimated by requiring x ~ H. 
Taking H to be constant, Eq. (l45p implies that t(x*) must solve 



In — = — In 

X* 4 



1 / n i \ 2 r i 3 / 2 

^(X,)J + ln ( v ^vV^- i ^^^J + 3[M(x*)J , (46) 

in which the Planck mass, Mp, has temporarily been restored to exhibit the relative 
magnitude of each term. The asymptotic regime of late times corresponds to \it ^> 1, 
and is approximately reached when N ~ H/fi e- folds have elapsed since the waterfall 
transition. In this region ln(/xt) can be neglected in comparison with (/zt) 3//2 . Under 
the same assumptions, (tMp) -1 can be neglected in the middle logarithm, after which 
the remaining term is of order ln/zMp 1 . This may be large if fi is much smaller 
than the Planck scale. On the left-hand side, however, In Hx* 1 is typically of order 
1.5 ln(A;*/ A;) + 0.51nMi/ _1 ~ 100. To avoid fatal problems with overproduction of 
black holes and topological defects, discussed by Garcfa-Bellido et al. [20J, the entire 
waterfall phase should complete within an e-fold. Since the scale for roll-down is set by 
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//, this requires \i > H. Hence /i cannot be arbitrarily small, and | ln^Mp 1 )! is likely 
to be much smaller than 100, and can be discarded in comparison with lni^x" 1 . The 
remaining terms imply that the transit time is given to a fair approximation by 

t(x*) = -(h*-) 2/3 - (47) 



A* V 2 X*J 

Note that the typical smallness of x*> which at first sight might cause one to imagine 
that any isocurvature fluctuation is tiny, has played an important role in ascertaining 
the rough validity of Eq. ( T4TI) . 

In virtue of our assumption that H is practically constant, Eq. ( )4"T|) can be converted 
to an estimate of the transit time in e-folds, N(x*)- The moments of N satisfy 

XT poo /o rr \ 2/3 

(N) = - P(x*,cr x )(-ln-) d **' ( 48 ) 
A* Jo V z a*/ 

<,v2 H7) I P{X "^\2^J «•*•■ < 49 > 
w = {f) [ PM ^S Ax " (50) 

where -P(x*> °x) ^ s ^ e distribution of initial values, x*- It is determined by the quantum 
mechanics of a massive scalar field evolving under Hubble damping. Therefore, we 
expect P to be close to Gaussian, characterized by its variance but with negligible 
higher moments. 

In analogy with the inverted simple harmonic oscillator considered in §D these 
expressions lead to the surprising result that the variance and skew of ( w depend only 
weakly on the initial variance, er^. Furthermore, even if the initial distribution P(x*, o~ x ), 
is exactly Gaussian, a skew of order unity is generated. This behaviour persists over an 
exponentially large range of <r^. We have evaluated Eqs. ( HHl) - (l50l) numerically, and used 
Eqs. f|T9l) - fr20|) to determine statistics of the observational quantity, ( w . Our results are 
depicted in Fig. [TJ To a good approximation, it is clear that (CwCw) = °"w ~ l{Hf fi) 2 , 
where 7 varies slowly within the range 0.01 < 7 < 0.1 over an exponentially large range 
of cr^. Similar behaviour occurs for the skewness, defined by 

Skew( Cw ) = (51) 
\Wsw/ 

and displayed in Fig. |2j Over a similar exponentially large range of a x it is clear 
that Skew(£ w ) is practically constant, taking values close to 1.5 for an Gaussian initial 
distribution of x*- 

We have verified these results using numerical Monte Carlo simulations, which 
give good agreement with the semi-analytic results given above. We have carried out 
simulations for a simple example of a potential which agrees with Eq. f )29|) in the 
transition region, but which is bounded from below and describes the descent into a 
stable reheating minimum at zero energy. The potential we use is 

A(0 2 -02) X 2- 



V = V n 



1 + 



4M 2 ; 4VWo 



(52) 
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~ Vo+ 2^/- \ VX 2 + 0(<p\ . . .) (53) 

where (p = <p — <Po- Near the hybrid transition, Eq. ( l53l) indicates that ( )52l) has the 
approximate form of Eq. (1291) . In using Eq. (153D we are free to end the simulation at 
any time after kinetic domination is reached. For convenience we choose the time at 
which the field is close to minimum for the first time. This occurs just after kinetic 
domination. We evaluate N as a function of the initial values 0* and x* by solving the 
full cosmological field equations. Initially is fixed so that the system is on the cusp 
of the hybrid instability. For the full potential, Eq. ( 152]) . this corresponds to choosing 
0* = <po- The initial value, x*i is drawn from a Gaussian distribution with variance 
a 2 ,. This process is repeated ~ 10 4 times, for which N is evaluated at a fixed value 
of H. This procedure accurately provides values for (N), (N 2 ) and (N 3 ) evaluated on 
a comoving hypersurface. Eqs. (Tl9l) - (l2"0j) are employed to determine the statistics of 
£. Physically, this procedure is equivalent to allowing an ensemble of trajectories, with 
initial conditions characterized by a distribution of values x — X*> to evolve under the 
influence of the hybrid potential. No approximations are made during the dynamical 
evolution. We include the effect of cosmological expansion on the evolution of and x-, 
and retain all contributions to the Hubble rate. 

We have performed this proceedure for a wide range of parameter values. The 
results are always in good agreement with our semi-analytic formulae. We find that 
the variance <j\ depends only weakly on the initial variance, cr^. More importantly, the 
skew is of order unity, irrespective of the value of a 2 . Our conclusion that ( w carries 
a strongly non-Gaussian distribution appears to be robust. In Figs. [TH2} we present 
results obtained using the full potential fl52l and the specific choices 

7] = -0.02, A = 10~ 5 , V = HT 11 , 0o = 1- (54) 

In this model H/fj, ~ 0.02. 

One may have harboured concerns the approximation of Eq. ( l32|) . in which <fi is 
taken to be constant, could be a good approximation near the top of the arete, but a 
poor description of the waterfall evolution near kinetic domination, where the dominant 
contribution to £ w is synthesized. The good agreement shown in Figs. [IH2] between 
our exact Monte Carlo simulations and the semi-analytic formulae shows that Eq. fl32l) 
provides a good approximation throughout the evolution of the waterfall field. 

3.3. Power spectrum and /nl 

According to the argument of § §3.11 - l3~2"l we have identified a new source of curvature 
perturbations in hybrid-like scenarios. However, we cannot yet conclude that this 
perturbation is present in typical hybrid models because the foregoing analysis depends 
on the validity of the separate universe picture. Our numerical experiments show 
that the dominant part of ( w is generated close to the era of kinetic domination. In 
the Introduction (JQ), we discussed the role of tachyonic preheating in quenching the 
dispersion of trajectories. This is associated with a failure of the separate universe 
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Figure 1. er^/(iJ//u) 2 plotted as a function of the RMS initial fluctuation a x in \- 
The green dots are computed numerically for a Gaussian ensemble of trajectories, using 
the potential (|52|) and parameter values (|54|) . The lines are semi-analytic estimates, 
obtained by numerically integrating (|48H50[) . The solid red line uses a Gaussian 
distribution of initial conditions of variance a x , while the dashed blue line employs 
a uniform distribution of varaince a x . The final variance depends weakly on the form 
of the initial distribution, and varies extremely slowly over many decades of the initial 
variance. 
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Figure 2. The skew of Cw plotted as a function of the RMS initial fluctuation a x in \- 
The green dots are computed numerically for a Gaussian ensemble of trajectories using 
the potential (|52|) and parameter values (|54|) . The lines are semi-analytic estimates, 
obtained by numerically integrating (|48H50[) . The solid red line uses a Gaussian 
distribution of initial conditions of variance a x , while the dashed blue line employs 
a uniform distribution of variance a x . The final skew depends weakly on the form 
of the initial distribution, and varies very little over many decades of variation in the 
initial variance. 
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picture. If tachyonic preheating completes before kinetic domination is attained on 
CMB scales, then the large-scale waterfall field acquires negligible dispersion. In this 
case the evolution of perturbations follows the traditional picture of Refs. [HJ ED] jjj] We 
study this possibility in detail in §HJ 

In the remainder of $3] we assume that the dynamics of some specific model are such 
that kinetic domination occurs before the onset of gradient instabilities associated with 
tachyonic preheating. In any such model ( w is uncorrelated with, and adds incoherently 
to, the curvature perturbation synthesized during the previous inflationary epoch. This 
is similar to the scenario studied by Boubekeur and Lyth [33], who considered the 
effect of an uncorrelated x 2 perturbation which added incoherently to a scale-invariant 
inflationary contribution. 

Does ( w unacceptably contaminate the primordial density perturbation in the 
hybrid scenario? The answer depends on the spectrum of ( w , and whether it makes 
a significant contribution on CMB scales. The spectrum of ( w is determined by the 
dependence of o\ on the moments which characterize the initial distribution of x*- The 
continuum power spectrum associated with £ w takes the form 

V (k) - - + V M^HL (KK) 

w{> dlnfc da 2 x din k + ^ 3 d(i m dink' 1 ' 

where /i m is the mth moment of x*i °\ is the initial variance of Xi an d V denotes 
the dimensionless power spectrum, defined for any fluctuation by the rule V(k) = 
k 3 P(k)/2ir 2 . Comparison of the Gaussian and uniform distributions in Fig. [T] (which 
have very different higher moments) gives reasonable evidence that da^/da^ is roughly 
independent of the /i m . There is no reason to believe that the dcr^/d/i m are small 
in comparison, but if the initial distribution of x is close to Gaussian then these 
moments can be neglected. As we argued above it seems reasonable to expect that 
P(x*) is approximately Gaussian, because it is determined by the fluctuations of very 
massive oscillators evolving under Hubble damping. Under these assumptions the power 
spectrum of ( w is approximately 



2 



V w (k) * P^V x (k). (56) 

x 

In view of the steep, blue spectrum associated with x* one might worry that Eq. (1561) 
implies V w is also very blue. This turns out not to be the case. A slowly- varying function 
such as cr^(cr^) can be well-approximated by a series in powers of ln(if/<r x ), of the form 



In 



(H 



A 1 



In 



( H 



x / J 



-i 



0(9") C,7! 



The parameter p must be strictly positive, because must be an increasing function 
of a 2 . Subsequent terms in the series are subleading. One can show analytically that 

|| A significant non-Gaussian perturbation is always generated on the horizon scale, which may or may 
not be competitive with the perturbations studied by Randall et al. [18] and Garcia-Bellido el al. [20] . 
We defer consideration of this effect until 35l 
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a 2 ^ takes precisely this form when the distribution for x* is uniform, with leading term 
equal to (2/3) 2 ^ 3 (H/n) 2 [\n(H/a x )]^ 2 ^ 3 . Focusing on the leading term, one obtains a 
power spectrum 

diner 2 ( n An) 

In 



K,w- dln<T ^- 4o) 



{- 



-p-1 

+ ■ ■ • . (5? 



din A; 2 

Eq. (|58|) shows that V w {k) is proportional to the tilt of the power spectrum of V x (k), 
rather then the input spectrum itself. Even if x nas an extremely steep power 
spectrum, provided its running is small, the power spectrum of £ w is approximately 
scale invariant. Had we retained subleading terms from (157j) . proportional to higher 
powers of [ln_ff/<7 x ] -1 , the same conclusion would have been obtained. Such terms may 
modify the normalization of V w but do not change its shape. 

Our next step is to estimate /nl- Our semi- analytic discussion and the numerical 
evidence shows that the skew of ( w is of order unity over a very wide range of a x , and 
is essentially independent of k. It follows that 

(CwCwCw) ~ (CwCw) 3/2 . (59) 

Assuming there are no other non-Gaussian contributions, /nl is given by the ratio 

f (CwCwCw) 'Pw(k) I (c.r\\ 

JNL - ^ 7T\2 ~ TM2' \ W 

V tota i{ky V t otA k ) 

where Ptotai is the total power spectrum of (, including the usual inflationary 
contributions as well as those from the waterfall process. The contributions from the 
waterfall need not dominate "Ptotai- Fig. [1] shows that when a x varies over ~ 30 orders 
of magnitude, o~ w varies only over a factor of roughly two. We will take a conservative 
point of view and estimate that the variance in ( w is given approximately by 



V w ~ (CwCw) ~ 7 y-J , with 7 ~ 10- 2 . (61) 

The dependence on H/fi follows from our semi-analytic arguments and is supported 
by our numerical evidence. The dimensionless factor 7 ~ 10~ 2 arises in the same way. 
From Fig. [T]it is evident that 7 varies from about 0.03 - 0.1 over many decades of o~ x , 
so Eq. flBTl) is conservative. In practice, 7 is slightly larger in all of the explicit examples 
we have studied. In what follows, we quote requirements on /nl in terms of constraints 
on H/ [l. At the end of this section we will relate these constraints to parameters of the 
inflationary model such as A, V, and e. 

In models where tachyonic preheating does not quench the production of ( w , the 
power spectrum and /nl put tight constraints on H/fi. Whatever the other predictions 
of a given model of hybrid inflation, the power spectrum of ( must match observation. 
Therefore its value is fixed at Ptotai ~ 3 x 10~ 9 . If V w makes a dominant contribution 
to 7-totai the model is experimentally ruled out, because a correct normalization of the 
primordial power spectrum implies /nl ~ ^totaf ~ 2 x 10 4 . This is in grave conflict 
with experiment [31], which requires |/nl| ^ 100. To avoid ( w dominating the power 
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spectrum, the estimate of Eq. (I6T1) shows that we must have V w < Ptotai, corresponding 
to (H/fi) < 5 x i0~ 4 . 

In the regime where ( w does not dominate Ptotai, a slightly stronger constraint on 
H/fi arises from considering / NL . Combining Eqs. (!60|) and (IBTj) . and neglecting the 
contribution of V w to Ptotai, we find 

/ NL ~10 14 f-V (62) 



V/i 

Since current experimental constraints require /nl ^ 10 2 , we obtain the constraint 

- < 10- 4 . (63) 
fi 

This is slightly stronger than the constraint from the power spectrum alone. If future 
constraints on /nl improve to |/nl| 1, an d we take a less conservative estimate that 
V w ~ lO^if/yu) 2 , the resulting constraint on H/fi would be an order of magnitude 
tighter. 

3.4- Constraints on hybrid model parameters 

The constraints on H/fi turn out to be surprisingly restrictive when written in terms 
of the fundamental inflationary model parameters A, V, and e. Using the Friedmann 
equation, the slow-roll condition for 0, and the definitions ( 129|) of A and ( 133]) of /i, we 
find 

1/3 

(64) 



H / V, 



end 



fi \6y/2M*\y/^ t 

Here, V^ n d and e en d are the potential and slow- roll parameter evaluated at the end of 
inflation, which corresponds to the beginning of the hybrid transition. (Recall also that 
e is defined using the derivative of V in the p direction only, just before the hybrid 
transition). In the following, we will find it convenient to parametrize V^j as 

Knd = Mf (65) 

where Mj is the inflationary mass scale, measured at the end of inflation. Since current 
cosmological scales have exited the horizon long before the end of inflation, there are no 
a priori constraints on M/ and e en d from experiment. 

We can get an idea of how strongly non-Gaussianity constrains the underlying 
hybrid inflationary model using a simple estimate, in which we make certain assumptions 
about the naturalness of the potential. On CMB scales the normalization of the scalar 
power spectrum constrains a combination of Mj and e, specifically 

24^ = *~ = 3Xl0 " 9 < 66) 

Because H and e are close to constant during a hybrid phase it would be unnatural 
for this ratio to change significantly by the end of inflation, although such behaviour 
could be introduced by a tuned or highly featured potential. Moreover, e does not 
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approach unity towards the end of inflation, because there is no graceful exit in a 
hybrid model. Therefore, we generically expect ( 166]) to hold, roughly, throughout the 
hybrid phase. (For example, it holds for the potential ( 152]) for reasonable choices of 
parameters.) Accordingly, we will take Mj ~ (3 x lCT^e^Mp. This immediately 
implies that /nl < 100 requires Mj < AlCT 8 Mp. Mj sets the scale of the inflationary 
potential, and even if A took a value close to the Planck mass this inequality would 
imply a rather low value for Mj. It is more reasonable to assume that A would take a 
value at a similar scale to Mi — otherwise we would be forced to introduce a hierarchy 
between these mass parameters. Taking A ~ Mj implies the constraint 

M 7 < KT 8 M P . (67) 

The /nl parameter scales linearly with Mj, so if future experiments constrain |/nl| ^ 1, 
this constraint tightens to Mj < 10~ 10 Mp. We conclude that unless the inflationary 
scale is very much lower than the Planck scale, unacceptable non-Gaussianity will be 
introduced. 

The constraints we have derived here should be viewed as guidelines for reasonable 
hybrid inflationary models. If one has a particular example in mind, one should first 
compute the dimensionless factor 7 appearing in (1ST]) — although in practice we find 
this parameter varies very little between different potentials. One must then apply the 
definition of /nl i n (|60|) to find the corresponding constraint on H/ fx. The formula ( 1641) 
constrains the underlying inflationary parameters V, A, and e. The resuting constraints 
may be summarized by 

.1.1715*) ~ IjFsJ {WmJ {-^) (l^J • (68) 

where we have assumed that H//i is small enough that the waterfall makes a 
subdominant contribution to the ( power spectrum. While there may be models for 
which Eq. (1681) is effectively unconstraining — especially hybrid models with a very low 
inflationary scale and a correspondingly tiny e — our analysis indicates that there are 
surprisingly strong constraints on a broad family of hybrid models. Moreover, these 
constraints can be expected to improve as future experimental results become available. 



4. The validity of the separate universe picture 

In this section we return to the question posed at the beginning of §3.31 namely 
whether the trajectories followed by adjacent cosmological patches can be prevented 
from dispersing sufficiently to generate an appreciable ( w . Felder et al. [25] and 
later Felder, Kofman & Linde [26] argued that a phenomenon known as tachyonic 
preheating typically occurs in hybrid models. Lattice simulations of the effect show a 
rapid destabilization of the background evolution owing to gradient instabilities formed 
as the universe reheats on small scales. This destabilization invalidates the separate 
universe assumption. To determine how much dispersion occurs, we must decide whether 
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CMB-scale regions of the universe reach kinetic domination before or after the onset of 
gradient instabilities. 

Which are the relevant timescales? To avoid overproduction of topological defects 
and black holes it is already known that the hybrid transition should complete in of 
order one e-fold [20]. This constraint is independent of £ w . Accordingly the shortest 
wavelength modes, which leave the cosmological horizon immediately prior to the 
tachyonic transition, must reach kinetic domination after < 1 e-fold. Their effective 
release point xl short ^ can be found using ( 1431) and ( 1451) with k = k±, giving approximately 



(short) = / CH3 m 

We recall that C is implicitly defined in Eq. (HTJ) . and is of the form C ~ a + biMj \i) 2 , 
where a and b are roughly unity. Together with ( H71) . Eq. ( 1691) allows us to estimate the 
number of e- foldings, N$, for the shortest wavelength modes to reach kinetic domination. 
It is 

Ar H f. 67r 3 M\ 2/3 

iV^-(ln w ] . (70) 

The analysis of Garcfa-Bellido et al. [20] shows that the logarithm cannot be larger 
than (if//i)~ 3//2 . We also define Ni to be the number of e-foldings required for CMB 
scales to reach kinetic domination. These scales left the horizon iV* e-foldings before the 
end of inflation, and hence the CMB scale is roughly e^* larger than the horizon scale. 
In typical models iV* is of order 50 - 75. The scaling of x* with k/k+ indicates that 

x (lon g)=e -3^/ 2x (s h o rt )_ (n) 

Hence, 

H ( 6vr 3 M\ 2/3 

In typical models, the evolving field rolls into its reheating minimum almost 
immediately after kinetic domination. When it does so, Felder et al. [25, 26J argued 
that a rapid growth in occupation numbers of low-lying Xk modes would efficiently 
drain energy from the rolling field, so that preheating would complete within a single 
oscillation. On this basis, we should expect the separate universe picture to provide an 
accurate description up to ~ few x Ns e-folds from the hybrid transition, but to fail at 
later times. A definitive determination of the range of validity of the separate universe 
picture will likely require a full lattice simulation of the waterfall phase. As a proxy for 
this unknown condition, we will assume that if Nl/Ns < 10 then there is a reasonable 
basis for belief that kinetic domination on CMB scales occurs sufficiently quickly for £ w 
to be synthesized. 

There are two regimes in which the estimate of N^/Ns is especially simple, 
depending on the relative size of the two parenthesized terms in (ITT!) . Although one 
might have expected that the exponential suppresion of CMB-scale waterfall fluctuations 
would imply Nl ^> Ns, we will find that in both regimes, there is good evidence that the 
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synthesis of ( w can be successfully completed. The parity between N L and N s improves 
when H/ /i is made small enough to satisfy observational constraints. 

The first regime applies when N* > \b.(6tt 3 M/HC). This regime is relevant when 
no large hierarchies exist between M, H, and //. In this case, 

N L - N s ~ — (3i\g 2/3 , ioiN, > ln(67r 3 M/i/C). (73) 

H 

Both Nl and Ns are much smaller than unity, since they are proportional to the very 
small parameter H/fi. In §3.31 we showed that consistency with constraints on /nl 
requires H/fi < 10~ 4 . For iV* ~ 50 — 75 this gives Nl — N s < 3 x 1CT 3 , showing that 
CMB-scale modes reach kinetic domination only a thousandth of an e-fold after short- 
scale modes. In this regime the ratio Nl/Ns is controlled by the hierarchy between N* 
and the logarithmic factor, yielding 



N L 
Ns 



1 + SN - 



\n(Qn 3 M/HC) 



2/3 

for JV, > ln(6vr 3 M/HC). (74) 



Eq. flTDl shows that the logarithm may be of order iV* without spoiling the requirement 
Ns < 1. In the typical models studied in $3] we have found that Nl/Ns is generally of 
order 5-7. This analytic estimate is confirmed by our detailed Monte Carlo simulations. 
Pending detailed lattice calculations, there seems no reason to believe that perturbation 
theory should fail long before kinetic domination occurs on CMB scales. We also 
note that, while inflationary initial conditions imply that arguments based on causality 
give limited information, it may happen that the transition from an inflationary to a 
radiation-dominated equation of state takes some time to propagate from short to long 
scales. If so, it may be possible to tolerate larger values of N L /N S than might seem 
reasonable on the basis of Refs. [251 EE] . 

Even if a typical member of a bundle of CMB-scale trajectories does not reach 
kinetic domination, a small subset of trajectories in the bundle can be expected to do 
so. In this case the final curvature perturbation is likely to depend in detail on the tails 
of the distribution of \*i an d cannot be calculated in perturbation theory. We leave the 
exact details of this regime for future work. 

The second regime corresponds to iV* <C ln(67r 3 M/HC), which requires large 
hierarchies between M, H, and fi. In this regime, Nl and Ns are very nearly equal, 
with 

H ( fivr 3 M\ 2/3 

N L ~ N s ~ - ( In -j^- J for N* < \n(Qn 3 M/HC). (75) 



It follows that Ns must be exceedingly small, since 

^^~ H Jm/HC) ^N,«M^M/HC). (76) 

Then, because N±/ ln(6ir 3 M / HC) -C 1 by assumption, the difference Nl — Ns is itself 
much smaller than iV^. In this regime there seems a very good basis for belief that the 
separate universe assumption is not invalidated. Note that, however large the logarithm 
becomes, the requirement Ns < 1 can be satisfied by choosing a correspondingly small 
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5. Conclusions 

We have argued that the transition ending hybrid inflation can generate a large 
contribution to the curvature perturbation. This contribution arises from fluctuations in 
the "waterfall" field, Xi which are created as modes exit the inflationary horizon. These 
perturbations are subsequently suppressed by the cosmological expansion. Although 
these fluctuations are very small when inflation ends, they are amplified by the tachyonic 
instability in \ and can provide a significant contribution to the curvature perturbation. 
To calculate the enhancement we use a gradient expansion in which we treat the 
moments of the initial distribution non-perturbatively. This method is similar to the 
method of "moment transport" which has recently been used to calculate the curvature 
perturbation in inflationary models [32]. Here, the procedure is simplified owing to 
our ability to calculate an approximate analytic solution. This technique predicts that 
the fluctuations generated by the waterfall are highly non-Gaussian, and can easily 
exceed observational constraints. If present in a given model, these effects lead to new 
constraints on viable hybrid models together with a unique observational signature. 

There are two phenomena which enable the tiny fluctuations in \ to generate 
significant effects. First, the variance in 5N following completion of the waterfall 
depends extremely weakly on the variance in \ before the waterfall begins. While 
surprising, this phenomenon is sufficiently general that it can be illustrated using an 
inverted simple harmonic oscillator. Similar behaviour is seen in the semi-analytic and 
numerical examples given above. Therefore, even very tiny fluctuations in \ can lead to 
large fluctuations in 5N. This means that the exponential suppression of ^-fluctuations 
during inflation does not guarantee that these fluctuations are harmless when hybrid 
inflation completes. In a single-field model, this kind of amplification would be forbidden 
because the curvature perturbation ( would be conserved on super-horizon scales. Since 
hybrid inflation necessarily involves multiple fields, the presence of isocurvature modes 
means that ( is no longer conserved. As we have argued above, the hybrid waterfall can 
cause it to grow significantly. 

Second, the spectrum of fluctuations in d~N, and hence the spectrum of the curvature 
perturbation (, is nearly scale-invariant. Naively, one might have expected that any 
process operating at the end of inflation could only affect modes of wavelength smaller 
than the horizon. This intuition is borne out for some processes: for example, the 
production of black holes or defects at the end of hybrid inflation results in a very blue 
power spectrum, which would be unimportant on cosmological scales. In the present 
case, it is true that \ has a strongly blue power spectrum immediately prior to the hybrid 
transition. This is because fluctuations in x are suppressed during inflation. However, 
the waterfall process ensures that the final fluctution in 5N is nearly independent of 
that in x- Therefore, the final curvature fluctuation is nearly the same, independently 
of scale, even though it is seeded by fluctuations which, on long wavelengths, are much 
smaller. Thus, the hybrid transition can reprocess the blue spectrum of x hrto a nearly 
scale-invariant one. The key difference in behaviour between large and small scales is 
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the transit time from the onset of the waterfall to kinetic domination. 

The analysis we have presented here depends crucially on the validity of the separate 
universe picture while CMB scales are being ejected from the arete. On small scales, 
it is known that the hybrid transition creates large (nonlinear) gradients in the scalar 
fields. This rapid growth in gradients is responsible for the "tachyonic preheating" 
effect studied by Felder et al. [251 126]. It is possible that short wavelength (horizon- 
size) perturbations reach this nonlinear regime very rapidly, before significant curvature 
perturbations are generated on longer wavelength (cosmological-size) scales. If so, this 
would prevent the dispersion of trajectories and lead to effects which are much weaker 
than those estimated here. An estimate of the onset of this gradient instability can be 
obtained by studying the time taken for different scales to reach the kinetic-dominated 
regime, at which time the curvature pertubation ( ceases to evolve. Our estimates 
indicate that long and short wavelength perturbations reach the kinetic-dominated 
regime at nearly the same time, within a tiny fraction of an e- folding from the onset 
of the hybrid transition. The hierarchy between the transit time for long and short 
wavelengths is typically a number of order unity. In this regime a subtle interplay exists 
between effects on short and long length scales, and the condition for validity of the 
separate universe picture is not known. We expect that a definitive resolution of these 
questions requires numerical simulations of the hybrid transition. We conclude that 
perturbation theory provides no compelling reason to believe that waterfall contributions 
to the curvature perturbation should be strongly suppressed. 

Even if tachyonic preheating prevents the generation of a significant curvature 
perturbation on CMB scales, it seems impossible to prevent the appearance of strong 
non-linearities on the shortest scales — those associated with the horizon scale at the 
transition. If these fluctuations are large, they will collapse to form black holes or 
topological defects. The analysis we have presented indicates that the formation of 
these non-perturbative objects should be taken to proceed from very non-Gaussian initial 
conditions. It is known that the final mass fraction contained in such collapsed objects 
can depend on the detailed profile of initial fluctuations [351 [36], [37] , which potentially 
provides another means to constrain models containing a waterfall transition. 
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Appendix A. Evolution of a massive scalar field in de Sitter space 

Here we give more details regarding the calculation of the power spectrum of % before 
the tachyonic transition. This is used to determine <J^(k), the variance of \ smoothed 
on a real space length scale of 2ir/k. 

We assume that x has the effective Lagrangian 



S 



-gd 4 x (A.l) 



and take m to be constant. We work in conformal time and assume that the spacetime 
is flat de Sitter space, and hence that 

^---^ <"> 

where H is the Hubble parameter. As is conventional, we take rj to increase along 
(—00, 0) as proper time increases and the universe expands. We will find it convenient 
to use a rescaled field u, defined by 

u{t,x) = {-H V ) X {t,x) (A.3) 

Substituting this definition into the action for x, and decomposing u into its Fourier 
modes, we find that each Fourier mode possesses the effective action 



S = M) 2 - J * 2 - - - «1 * (A.4) 



where ' = d/drj. The rescaled field u has a canonically normalized kinetic term in this 
action. This makes it simple to impose the usual Minkowski space boundary conditions 
well inside the cosmological horizon. The subtlety is that well outside the horizon the 
mode functions are very different from their Minkowskian form. 
The equation of motion resulting from ( 1A.4I) is 



with solution 

uM = c iy pk^HS-\-kri) + c 2 ^k^H^\~kn) (A.6) 
where H^' 2 ^ are Hankel functions and 



For applications to hybrid inflation, we require that the waterfall field is very massive, 
and in particular that M ^> H. Hence the index of the Hankel function is purely 
imaginary. We assume that we work in this limit throughout this derivation and hence 
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take v = i(M/H). We now canonically normalize the mode functions by taking the 
limit — krj — >• and using the small argument expansion^] 



H^' 2 >(x)^\ — exp I ±ix =F — — =F — ] x > \v\ > 1 (A.8) 

V 7TX \ 2 4 / 



To canonically normalize the fields, we should have the positive frequency mode 
asymptote to Uk(x) — > (2k)~ l l 2 e~ lx where x = —krj. Keeping in mind that /i is 
purely imaginary, this implies that the correctly normalized positive frequency mode 
is described by 

e-M/2 ft 



To find the super- horizon limit, we take x <C 1 and use 

x < 1 (A.IO) 



7T 



' (£)" i» (| 



Starting with the Euler reflection formula 

r(-l/)]» = r^— (All) 

sm(7rz/) 

and using the fact that T(z/)* = T(z/*), the definition of x m terms of u, and the 
assumption that v 3> 1, the super horizon limit of x becomes 



Xfc ^ il^^ k ^ /2)+miH)+i6 ^ >>1 ' M /^»l- (A.12) 
where # is an irrelevant constant phase angle. Using this expression, one can easily 
derive the correlators (j4T)l) of Xk on super-horizon scales. 
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